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1. Introduction 

Brownian motors are small scale engines that operate on the basis of nonequilibrium 
fluctuations. Most of the models that have been discussed previously in the literature 
refer to the transfer of chemical or potential energy into work [1] by Brownian entities 
that operate in a externally imposed asymmetric environment (or by spontaneous 
symmetry breaking [2]). In this 100 percent efficiency can in principle be 

achieved [3], in accordance with the basic laws of thermodynamics. Recently, a thermal 
Brownian motor was introduced [4], for which an exact microscopic analysis turns out 
to be possible [5]. The model can be viewed as a simplification of the Feynman- 
Smoluchowski ratchet [6, 7]. This construction differs in two essential aspects from 
most other Brownian motors. First, the source of energy is the presence of different 
heat baths. Thermodynamics thus predicts that the limiting efficiency is that of the 
Carnot engine. Second, the asymmetry is an intrinsic property of the motor, not of 
the environment or of the imposed constraint. The purpose of this paper is to review 
some further properties of this construction, and to discuss in particular specific issues 
including the relation with macroscopic response and the efficiency at maximum power. 



2. Model and basic results 



We first briefly review the basic construction of the thermal Brownian motor and the 
main results. The motor consists of a set of convex units i (i = 1,...,N), that move 
as a rigidly linked entity with a single degree of freedom, along a specified direction x. 
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J dV W(V\V')P(V',t) - W(V'\V)P(V,t) . (2) 



Each unit resides in a separate, infinitely large compartment that is filled with an ideal 
gas at equilibrium. These compartments play a role akin to the thermal reservoirs in 
a Carnot construction. For simplicity, we consider the case of spatial dimension d = 2 
with coordinates (x, y) . In Fig. 1 we show four different realizations of the motor. The 
gas particles are point particles of mass m. Their density, temperature and Maxwellian 
velocity distribution in compartment % are denoted by p iy Ti and 0j respectively, with: 

= 2^M ^ { 2k B T t J' (1) 
We consider the limit in which the mean free path of the gas particles is much larger 
than the linear dimension of the motor unit. In this case, there are no pre-collisional 
correlations between the speed of the motor and that of the impinging particles. The 
probability density P(V, t) for the speed V = (V, 0) therefore obeys a Boltzmann-Master 
equation: 

dP(V, t) 

di 

VF(T/|y) = * s the transition probability per unit time for the motor 

to change speed from V to V due to the collisions with the gas particles in the 
compartments % — 1, . . . , N. Its explicit form, following from the collision rules and 
basic arguments from kinetic theory, is given by (see [5] for more details): 

/*27r /*+oo /*+oo 

Wiiy\V') = / d9S i F i (9) / dv' x / dv' y pM<,v'y)(V' - if) ■ e ± 

Jo J —oo J— oo 

x e \{v> - ^ ■ ex] s\v - v - ^r'ln K - < cot d - v ')]- ( 3 ) 

L J l 1 + j± sin (J y 

Here O is the Heaviside step function. The 5 Dirac function selects appropriate 
particle velocities yielding the post-collisional speed V for the motor in function of 
the precollisional speeds V and v' x , v' y of motor and gas particle respectively. The shape 
of the motor unit i is characterized by a normalized shape function Fi(0), such that 
Fi(9)d9 is the fraction of the outer surface of the unit that has a tangent between 9 
and 9 + d9, cf. Fig. 2. The angle 9 is measured counterclockwise from the horizontal x 
direction. The total perimeter will be denoted by Si. 

Solving the Boltzmann-Master equation in general is out of the question, even at 
the steady state. To make progress, a perturbative solution is necessary. Since we 
expect that the rectification disappears in the limit of a macroscopic motor, a natural 
expansion parameter is the ratio of the mass m of the gas particle over the mass M of 
the object. More precisely, we have used e = ^Jm/M as the expansion parameter. 

Introducing the transition probability Wi(V';r) = Wi(V\V) with jump amplitude 
r = V — V, the Boltzmann equation can be rewritten as 

dP ^^ = j ' W{V -r;r)P(V - r,t)dr - P(V,t) j W(V;-r)dr. (4) 

The Taylor expansion of the r.h.s. of this equation with respect to the jump amplitude 
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Figure 1. Four different realizations of the thermal Brownian motor, each consisting 
of two rigidly linked motor units, (a) will be referred to as Triangula and is composed 
out of two identical triangles, while model (b) is called Triangulita and has a triangle 
and a bar. Constructions (c) and (d) consist respectively out of two identical triangles 
pointing in the opposite direction and two triangles with a different apex angle. 



e + de 




Figure 2. A closed and convex object with perimeter S. The length of the surface 
with an orientation between 9 and 9 + d9 is SF(9)d9, defining the form factor F(9). 



leads to the equivalent Kramers-Moyal expansion: 



dP(v,t) ^(-i) n fd\ n r , rrs „ /rr ^ 

n=l x ' 

with jump moments a n (V) = Yli a n,i(y)i defined as: 

a n ,iiV) = I r n Wi(y-r)dr. (6) 



Thermal Brownian motor 



4 



One can obtain from the Kramers- Moyal expansion the set of coupled equations 
for the moments of the velocity distribution (V n ) = J V n P(V,t). Of interest for the 
subsequent discussion, we reproduce the exact results for the first two moment equations: 



d t (V) =( ai (V)) 

d t {V 2 ) = 2(Va l (y)) + {a 2 (V)), 



(7) 
(8) 



with 



a lti (V) = —piSi PdBF0) -— 1 , 
J M + m si 



sin 2 9 

2k B Tim m 2 . 2 . 2 

exp| — — — —V sin 9\ V sin 



7T 



2k B T,i 



+ 1 + Erf 



m 



r>27T 



a 2) i(^) = / d9F i {9) 

Jo 



2k B Ti 
sin 2 6 1 



l^sin^ 



(k B Ti sin 6 + mV 2 sin 3 6) 



(9) 



2k B T(m m 
exp[ 



7T 



(M + msin 2 ^) 2 

\/ 2 sin 2 6} (2k B Ti + m\/ 2 sin 2 6) + 



2k B Ti 



+ 1 + Erf 



2k B T t 



Vsin9 



(3k B T im V sin + m 2 V 3 sin 3 0) 



(10) 



3. Systematic speed: general results 



The equation for the first moment is coupled to the higher moments, a familiar feature in 
statistical mechanics. However, a perturbative solution can be obtained by considering 
the Taylor expansion in e. To lowest order in e the first moment equation reduces to a 
closed equation featuring a linear relaxation law, namely: 

Md t (V) = - 7 (V), (11) 

where 7 = X^7j an d 1% the linear friction coefficient, due to the section of the motor 
sitting in gas mixture i: 

'k B T im 



7i = ApiSi 



2tt 



d9Fi(6) sin 2 9. 



(12) 



Note that at this level of the perturbation, the steady state speed of the motor is zero: 
even though the motor may be constituted of asymmetric parts, no rectification takes 
place at the level of linear response. 

At the next order in the perturbation, the first moment is coupled to the second 
one, which needs to be solved at lowest order. By doing so, one finds at the steady 
state: 



M(V 2 ) = k B T eS , 



(13) 
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Figure 3. (a) Top: The response function ai^(V), see Eq. (9), for different geometric 
objects used in de constructions of Fig. 1: a triangle with 6$ = 7r/36, a bar of length 
L and a disk with radius R. 

(b) Bottom: The function a\.i{V) for a triangle with opening angle respectively equal 
to 7r, 7r/2 and 7r/20. The following parameter values were used: mass ratio M/m = 20, 
density pi — 0.0022, temperature Ti = 1, size L = R = 1 and fee = 1 by choice of 
units. 



with T e fj the effective temperature given by: 

T eff = (14) 

By inserting this result in the equation for the average velocity, one finds for the lowest 
order term of the steady state velocity V : 

VmV 8M ^S^f'-MFMvAO) ' 

This speed is equal to the expansion parameter times the thermal speed of the motor and 
further multiplied by a factor that depends on the geometric properties of the object. 
Note that the Brownian motor ceases to function in the absence of a temperature 
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difference (when Tj = T c g for all i) and in the macroscopic limit M 



oo 



since 



Vq ~ \/M). Note also that the speed is scale- independent, i.e., independent of the 
actual size of the motor units: Vq is invariant under the rescaling Sj to CS{. In Table 1 
we have collected the explicit result for V for the four different thermal motors, depicted 
in Fig. 1. 

It is tempting to interpret the direction of the systematic speed on the basis of 
intuition. One can study the response of the motor by the application of an external 
force. Due to the asymmetry, the response curve (speed versus force) will be asymmetric. 
In particular, it will require less force to move a triangle in the "easy" direction of the 
sharp angle. Consequently, one expects that fluctuations leading to an "extra force" 
in the easy direction would be more successful in moving the object than the ones 
in the other direction. A similar argument can be applied to the discussion of the 
electric current in a diode, which also has an asymmetric response of current versus 
field. One expects that fluctuations in the electric field will produce a larger current 
in one direction than in the other one. However, as is well known [8], these arguments 
are clearly wrong. The first observation is that there is no systematic speed (or electric 
current) when the system is at equilibrium. In fact, the effect of possible asymmetries in 
a system are completely wiped out in equilibrium, as expressed explicitly in the principle 
of detailed balance discovered by Onsager [9]: at equilibrium any process and its time- 
reverse occur equally frequently. Second, equilibrium is typically a so-called point of flux 
reversal [7]. For example, in the case of Triangulita, the flux is in the "easy" direction 
of the sharp angle, only in the case when the temperature is lower in the reservoir 
containing the triangle. Otherwise the speed is in the opposite direction. Furthermore 
the speed again becomes zero when the temperature is zero in the reservoir containing 
the flat paddle. Finally, the above mentioned asymmetry in the response appears in 
the nonlinear regime. In this regime, the dynamics of the various moments are however 
coupled. In particular, the equation for the average speed is not closed, but depends 
on the entire probability distribution. The nonlinear response to an external force is 
therefore not linked in an obvious way to the intrinsic nonlinear dynamics generated by 
nonequilibrium fluctuations. 

Full information on the equation of the first moment is contained in the nonlinear 
function ai^{V). It is reproduced in figures 3(a) and 3(b) for the motor units used in the 
constructions of Fig. 1. However, the direction of motion, to lowest order in e, contains 
much less information. Indeed, from Eq. (15), one finds for a motor with 2 units: 



Hence the direction of the average speed does not depend on the densities of the gases 
or the size of the motor units. The quantity — 1 < (sin 3 9) x = J Q 27r d9Fi(9) sin 3 9 < 1 is a 




(16) 
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Stationary velocity V (order e) 


Triangula 


1(a) 


V2 :t Bm ( 1 sin 9 ) Pi^(T 1 -T 2 )(V7r-V7¥) 


Triangulita 


1(b) 


^:if m (l sin 2 g ) r ggl^glzZk) 


Triangle - triangle (><) 


1(c) 


<j2nk B m{i a \ piP2(T 1 -T 2 )(y/J\+y/7 2 ~) 


Triangle - triangle {6\ ^ 9 2 ) 


1(d) 


V27rfe B m pip2(l+sin6»i)(l+sin6l2)(Ti-T2)[v / Tr(l-sin 2 )— V^l-sinfli)] 


4M [piv/TTCl+sineO+pav^Cl+sin^)] 2 



Table 1. Analytic result for the lowest order contribution to (V) for the different 
constructions depicted in Fig. 1. 



measure for the sharpness of the unit: 

{— > 1 infinitely sharp to the left 
= symmetric upon reflection about y-axis (17) 
— > — 1 infinitely sharp to the right. 

On the other hand the quantity < (sin 2 6) < 1 is proportional to the linear friction 
coefficient, cf. Eq. (12). For a bar of zero width, it has the minimal value (sin 2 6) = 
when the bar is oriented horizontally, while it is maximal ((sin 2 #) = 1) for a vertical 
bar. 

Based on Eq. (16), one can distinguish four different classes of motors, depending 
on how the direction of speed is affected upon interchanging the temperatures T\ and T 2 . 
The examples of Fig. 1 each represent one such class. Fig. 4 illustrates their behavior 
upon switching 7\ and T 2 . 

4. Systematic speed: special cases 

To gain further insight on the direction of the systematic speed, it is instructive to 
calculate its explicit value for a number of special cases, see also Table 1. 

4-1. Identical units and nonlinear response 

For a construction with identical motor units, see for example Triangula (Fig. 1(a)) 
with two triangles pointing in the same direction, Eq. (16) reduces to sgn(Vo) = 
sgn(— (sin 3 9)). One concludes that the direction of Triangula's motion is independent of 
the temperature gradient. This behavior originates from the permutational symmetry 
of the identical motor units, implying that the speed must be invariant under the 
interchange of Ti,pi with T 2 ,p 2 - Note that the average velocity is pointing in the 
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Figure 4. Plot of the average speed (see Table 1) for the four different models depicted 
in Fig. 1 as a function of the applied temperature difference. The following parameter 
values are used: M/m = 20, pi = p 2 = 0.0022, O = tt/36, 6 x = tt/1000, 9 2 = tt/3 and 
T 1 +T 2 = 2 fixed. 
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direction of the sharp angle of the motor units. It turns out that this feature is general: 
for a motor with two identical units, Eq. (15) simplifies to 



V2^k^ (T 1 -T 2 )(V7\-Vn) 



sin 



3 



v ° ~ [P1V ^ +WV ^] 2 (^oy (18) 

We conclude that the motion is always in the same direction, the "easy" direction 
determined by the sign of (sin 3 #), cf. Eq. (17). 

We again stress that while the above direction of motion appears to be plausible, 
in agreement with our intuition, this specific feature is exceptional and a result of the 
permutation symmetry of the units. In fact, the latter symmetry has a dramatic effect 
on the response of the motor to the temperature gradient: as is clear from Eq. (18), the 
speed has a parabolic profile in function of the applied temperature difference T\ — T 2 . 
In particular there is no linear response regime, with a speed proportional to 7\ — T 2 ! 
It is due to this unusual, symmetry induced property that equilibrium is not a point of 
flux reversal. 

4-2. Identical units, with opposite orientation 

Another interesting class of motors are those with identical units, but pointing in 
opposite directions, cf. Fig. 1(d). For these motors is F\{9) = F2(tt—0) and consequently 

(sm n 6) ={ {sinn0) > n6Ven (19) 
{S h \ -(sm n 6) 2 nodd. [ ^ 

One finds from Eq. (15) that the average velocity is given by 
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We conclude that the equilibrium state 7\ = T 2 is the (only) point of flux reversal. In 
the particular case of equal density (pi = p 2 ), we have a symmetric Brownian motor, 
i.e. the absolute value of its speed is invariant upon inverting the temperature gradient, 
cf. [10]. 



4-3. Different units 

When the shape of the two motor units is different, the situation can be quite more 
involved. In particular, there may be multiple points of flux reversal in function of 
the applied temperature difference. This reinforces the statement that the direction 
of motion is difficult to predict. As an explicit example, we show in Fig. (4) that a 
construction with two different triangles (Fig. 1(d)) has a second reversal point for 
T\ > T 2 . 



4. 4- Multiple reservoirs 



To illustrate the effect of multiple reservoirs on the amplitude and direction of the speed, 
we consider the case with three units, namely two triangles (with different orientation 
and apex angle) and a bar, see Fig. 5. The stationary velocity reads 

V / 27rm&B 1 



V = 



4M [ 2pi L x v^T + p 2 L 2 ( 1 + sin 6 2 ) VT 2 + p 3 L 3 ( 1 + sin 9 3 ) : 



x 



p 2 (l - sin 2 9 2 ) 2p 1 L 1 ^T 1 (T 2 - 7\) + p 3 L 3 (l 



sin 1 



+ p 3 (sin 2 9 3 - 1) 2 Pl L 1 ^/T 1 (T 3 - T\) + p 2 L 2 (l + sin# 2 )- 



VT 3 (T 2 -T 3 ) 
^(T 3 - T 2 )] } . 



(21) 



The speed is of the same order of magnitude as before. The net motion disappears at 
equilibrium (Ti = T 2 = T 3 ), when the asymmetry vanishes (9 2 = 3 = 0, both triangles 
are flat), or when the effect of the two units cancels each other (for example 9 2 = 6 3 , 
T 2 = T 3 and p 2 = p 3 , i.e., both have the same opening angle in a reservoir at the same 
temperature and with equal density). 

One may wonder what is the effect of the addition of a reservoir that is at 
equilibrium with a motor consisting of two units (at lowest order in e). More 
precisely, consider Ti to be equal to the effective temperature of reservoirs 2 and 3, 
Ti = ( 72 T 2 + 7 3 T 3 )/(72 + 7s)- We find from Eq. (15): 

V2^h~k^ p 2 S 2 (T 2 - TO (sin 3 G) 2 + p 3 S 3 (T 3 - T 1 )(sin 3 6) 3 



V 



(22) 



4M EtiP^VTKsm 2 ^ 
We conclude that the presence of motor unit 1 slows down the systematic speed. Its 
single effect is the addition of an extra friction term, cf. the denominator in the above 
formula. 
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Figure 5. Construction of a thermal Brownian motor with 3 reservoirs. The three 
motor units are rigidly linked and move as one single entity. 

5. Efficiency 

As pointed out in the introduction, the efficiency of a thermal motor, like the one 
discussed here, has to be compared with that of a Carnot engine. The fact that we can 
extract work by adding a load to the motor implies that there must be a concomitant 
transfer of energy from the hot to the cold reservoir. A fundamental difference with the 
Carnot construction is however that the motor is at all times simultaneously in contact 
with both the thermal reservoirs (we consider for simplicity the case of two reservoirs). 
As a result, we expect that the object will acquire a temperature intermediate between 
T\ and T 2 . This observation is confirmed by the calculations mentioned above, where 
it is found that, to lowest order in e, the velocity distribution of the motor is Gaussian 
at an effective temperature T e g given by Eq. (14). At this order of perturbation, one 
can easily evaluate the corresponding heat flux between the reservoirs and show that it 
obeys a Fourier law, see below. Furthermore, it has been a matter of debate whether 
this irreversible entropy producing process is reducing (dramatically) the efficiency of the 
motor as a thermal engine. We will not discuss this matter further, see [11], [12], [13] and 
[4] for more details . We rather focus on a more practical definition of efficiency, namely 
the one considered in finite time thermodynamics. Indeed, in practical applications, 
Carnot efficiency is a less relevant quantity, since the central issue is not work but power. 
Hence, we will evaluate the efficiency of the Brownian motor at maximum power. To 
do so, we have to repeat the above analysis but including an explicit calculation of the 
heat flux and applying a (constant) load force F to the motor. 
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To explicitly include the energy transfer between the heat reservoirs, we add in the 
Boltzmann-Master equation, an additional variable E\ corresponding to the energy 
(change) in reservoir 1. Adding the direct acceleration produced by the external force, we 
obtain the following equation for the probability density P(V, Ei,t) (prime is denoting 
precollisional values) : 

dtP^E^t) = j dV'W^VlV^PiV^E'^t)- j dV'W 1 (V'\V)P(V,E 1 ,t) 

+ J dV'W^V^PiV'^E^t) - J dV'W 2 {V'\V)P(V,E 1 ,t) 



(23) 



with W(V \V) = Wi(V'\V) the transition probability per unit time for the motor to 
go from speed V to V, given by Eq. (3). Conservation of energy in reservoir 1 implies 



-MV 2 + E x = -MV' 2 + E\ . 
2 2 1 



(24) 



For the transition probability Wi(V;r) = Wj(V'|V) with jump amplitude r = V — V , 
the conservation of energy permits us to rewrite the above equation as: 

d t P(V, E u t) = J drW 1 (V - r; r)P(V - r, E x + ^Mr(2V -r),t) 
- J drW 1 {V;r)P(V,E 1 ,t) 

+ J drW 2 (V -r;r)P(V - r, E u t) - J drW 2 (V;r)P(V, E u t) 



(25) 



Equivalently, we have: 
d t P(V,E 1 ,t)= [ dr 



8 m 2V +r) 



e 9v e 2 



8E 1 _ I 



+ 



e r av - i 



W 1 (V;r)P(V,E 1 ,t) 



d 



W 2 (V;r)P(V,E 1 ,t) - — l-P(V,E u t) ) . 



(26) 



Since (£1) = f dV f dEiEiP(V, Ei, t), the energy exchange between the reservoirs obeys 
the following equation: 



dtiE, 



dr I dV j dE±E 1 
(-r) n d r - 



™ (- r )n gn °^ (y(2V + r)) m d 



n\ dV n 

n=0 m=0 



ml 



dE^ 



+ 



n=0 



n\ dV n 



W^V-r) 
(27) 



The only non- vanishing contribution comes from the term n = 0, m = 1, whence 

f Mr 

dtiEj = - { I dr—(2V + rjW^V; r)). (28) 
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Using the jump moments a n ^{V) denned by Eq. (6), we can write for the heat flux from 
reservoir 1 to reservoir 2: 

M 

Q^ 2 = -ME,} = M(Va hl (V)} + —(a 2A (V)). (29) 

Again, the functions a^i and 02,1, cf. Eqs. (9)-(10), can be expanded in terms of 
e = ^Jm/M to find the corresponding expansion for the heat flux in the limit of large 
mass M. To lowest order, one finds that the energy flux obeys a Fourier law 

Q^ 2 = «m -T 2 ) + , (30) 



with 



kp 7i72 /Qi \ 

M 7l + 72 



5.2. Power 



The time evolution for the average speed obtained from the Boltzmann equation (23), 
d t (V) = ( ai (V)) + ^, (32) 

differs from Eq. (7) by the acceleration resulting from the addition of the force F. 
To lowest order in e, the steady state solution of Eq. (8) for the second moment is 
unaffected by F (i.e., M(V 2 ) = hsT^) and we find that 

MV)) = -^((V)-V ). (33) 

Here is 7 the friction coefficient, cf. Eq. (12), and V the average steady state velocity 
in absence of an external force, see Eq. (15). The steady-state solution of Eq. (32) thus 
reads: 

(V) = V + -. (34) 

7 

Note that F stop = — 7V0 is the stopping force (at this order of perturbation). Obviously, 
in order to extract work from the Brownian motor we need to apply a force smaller than 
F st0 p. In particular one has that F ~ 1/M. 

One defines the efficiency of the motor as the delivered work per unit time over the 
heat exchange per unit time (for simplicity, we assume from here on that Ti > T 2 ) 

„=^ = Z£W (35) 

Ql->2 Ql^2 

Since to lowest order of the perturbation in e the heat flux is not influenced by an 
external constant force F, we can maximize the efficiency by maximizing the power 
—F(V). One easily verifies that the latter is maximal when we apply a force equal to 
half the stopping force, F = — 7V0/2. The resulting efficiency of the thermal engine at 
maximum power is then given by: 

Ak B 7i7 2 Ti - T 2 
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Efficiency 77 at maximum power 
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tt ran ™ fl ^2 (v'TT-v^) 3 
32 M ^ fmC7 ; y / Tl T 2 (^/T 1 +^/T 2 ) 

w m (l-sin 2 o ) 2 /tT Ti-T 2 

16 M l+sin(9„ Y T 2 (2VTT+(l+sin6» )v / 72) 2 

t m/i • /) \2 Ti-T 2 
32 Ml 1 Smt/ 0j V ^ T7¥ 

Trmn , ^flUl 1 ™„/3^i-T 2 rv^(l-sin0 2 )-yi¥(l-sm6>i)l 2 
32Ml i + Sm ^l i + Sln ^ v ^ [^(l+BinffO+^l+smOa)] 



Table 2. Analytic result for the lowest order contribution to the efficiency rj at 
maximum power, in the situation of equal densities in the reservoirs (pi = P2, T\ > T2). 



One can immediately deduce from Eqs. (12) and (15), that the efficiency r] is very 
low, due to a prefactor m/M: for small m/M the heat conductivity is the dominant 
factor and both the rectified motion and resulting work are a second order effect. This 
observation is of course restricted to our discussion of the small e limit. In fact much 
higher efficiencies can be in principle be achieved when one moves away from this limit 
[14]. 

Since the expression for the efficiency of an arbitrary motor, Eq. (36), is rather 
complex, we have reproduced the efficiency for the models of Fig. 1 in Table 2. In 
particular, rj simplifies significantly for identical conditions in both reservoirs (equal 
densities p 1 = p 2 and equal objects Fi(9) = F 2 (9) and Si = S 2 ): 

tt m / (sin 3 0) \ 2 (VTi-VT 2 )* 
1 32M\(sm 2 9)J ^T{T 2 {^Ti + VT 2 )' 1 ) 
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